Taming chaos by impurities in two-dimensional oscillator arrays 
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The effect of impurities in a two-dimensional lattice of coupled nonlinear chaotic oscillators and 
their ability to control the dynamical behavior of the system are studied. We show that a single 
impurity can produce synchronized spatio-temporal patterns, even though all oscillators and the 
impurity are chaotic when uncoupled. When a small number of impurities is arranged in a way, that 
the lattice is divided into two disjoint parts, synchronization is enforced even for small coupling. The 
synchronization is not affected as the size of the lattice increases, although the impurity concentration 
tends to zero. 



I. INTRODUCTION 

Coupled arrays of oscillators are studied extensively 
in many fields of science because of their prevalence in 
nature. They are used as models for coupled arrays 
of neurons [Q, chemical reactions coupled lasers ||] 
or Josephson junctions Q|, charge-density- wave conduc- 
tors crystal dislocations in metals |6|, and proton con- 
ductivity in hydrogen-bonded chains |7]. Various models 
and coupling schemes have been proposed and analyzed 
previously A particular class are arrays of coupled 
oscillators, which exhibit chaotic motion when uncou- 
pled. This class includes the forced Frenkel-Kontorova 
model 1^ , which finds a straightforward physical realiza- 
tion in an array of diffusively coupled Josephson junc- 
tions [ p^pr] , in which the applied current of each junc- 
tion is modulated by a common frequency. The pos- 
sibility to obtain synchronized motion in such systems 
has been investigated recently by Braiman et al. for the 
case of one- (ID) and two-dimensional (2D) chaotic ar- 
rays of forced damped nonlinear pendula [ [l2[ and cou- 
pled Josephson junctions p^ . They observed the emer- 
gence of complex but frequency-locked spatio-temporal 
patterns, in which the chaotic behavior was completely 
suppressed, when a certain amount of disorder had been 
introduced by randomizing the lengths of the pendula. 
In Ref. ijl^ the same phenomenon has been investigated 
from a completely different point of view: It was shown 
for ID arrays of coupled chaotic pendula, that introduc- 
ing a single impurity at a particular site is sufficient to 
lead to complete synchronization. 

In this paper we study 2D arrays of coupled chaotic 
pendula. We ask for the minimal coupling as well as 
the influence of concentration and arrangement of impu- 
rities, needed to observe spatio-temporal patterns. Al- 
though for geometrical reasons one might expect that a 
single impurity cannot play the role it plays in ID ar- 
rays, we find that it is able to tame the chaotic behavior 
of an arbitrarily large 2D array, provided that the cou- 
pling is strong enough. A single impurity can produce 
synchronized spatio-temporal patterns, even though all 
oscillators and the impurity are chaotic when uncoupled. 



In order to observe spatio-temporal patterns for smaller 
couplings, the geometrical arrangement of the impurities 
is shown to play a crucial role. Specifically we show that 
an impurity configuration that divides the lattice into at 
least two disjoint parts is most appropriate for the cre- 
ation of synchronized solutions. Such configurations will 
always produce patterns that locally can be identified as 
lines, provided the coupling is above a threshold value. 
The resulting spatio-temporal patterns are stable with 
respect to an increase of the system size, indicating that 
an impurity concentration that tends to zero may suffice 
to lead to synchronization. 



II. MODEL 

We will focus our analysis on the model examined in 
Refs. 
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where n,m = 1,2. ..TV. Thus, there is a damped, driven 
pendulum with unity mass and length Z„„j on each site 
(n, m) of the lattice, subject to an ac and a dc torque. 
The parameters used are the gravitational acceleration 
g = 1, the dc torque r' = 0.7155, the ac torque r = 0.4, 
the angular frequency u — 0.25, and the damping 7 = 
0.75. Neighbouring pendula are coupled via a discrete 
Laplacian, where k denotes the coupling strength. We 
have chosen free boundary conditions, i.e. 9a,m ~ ^'i.m, 

^Af,m = dN+l,m,On,0 — &n,l, Qn,N — Qn,N+l and USCd 

a fourth order Runge-Kutta routine with a time step 
dt — 0.01 to numerically integrate Eq. (|l|). We care- 
fully checked that decreasing the time step to dt = 0.001 
did not alter our results. 

A very convenient measure that allows a quick visual- 
ization of the average global spatio-temporal behavior of 
the lattice, is the average velocity 
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at multiple times of the driving period T = l/ui. Us- 
ing this quantity to obtain a bifurcation diagram 
not only can ascertain if chaotic or periodic behavior 
is obtained, but in addition helps to identify the max- 
imum period of a pattern: Computing a[t) at each pe- 
riod t = jT of the driving, will lead to a periodic se- 
quence (Ti , . . . , CTp , iJi , . . . , when transients have died out 
and a spatio-temporal pattern of periodicity p ('Pp at- 
tractor') has emerged. In practice we inspect the last 20 
values of a-{jT) with j = 1, ...170, so that transients have 
died out. Thus, P20 attractors or attractors of larger pe- 
riodicity are not recognized as such but rather appear as 
chaotic attractors. 

Applying this strategy to an isolated pendulum, a bi- 
furcation analysis with respect to the pendulum length I 
was performed in Ref. jl^ . This approach revealed that 
each isolated pendulum is chaotic for values I = 1 ±0.002 
and that three more chaotic windows exist: Two nar- 
row ones at I w 0.84, and / w 0.52, and a broad one 
for I < 0.35. Thus we know whether a chosen pendu- 
lum with length I is chaotic or not. Furthermore, for the 
above chosen parameter values it is known, that pendula 



with L 



> 1 and 



< 1 show libration and rotation. 



respectively, apart from the windows, where chaotic mo- 
tion appears. 



III. RESULTS AND DISCUSSION 

The simplest configuration of a 2D lattice described 
by Eq. (^, is that of a single impurity in a sea of iden- 
tical chaotic pendula with length = 1. We fixed the 
lattice size to be 50 x 50, with the the impurity located 
at site (25, 25) and made a bifurcation analysis with re- 
spect to the impurity length ^imp for various values k 
of the coupling. In Fig. |l| the obtained bifurcation dia- 
grams for fc = 1, 2, 3, 5 are plotted versus hmp- For con- 
venience we normalized a to take on values in the unit 
interval [0,1]. Similarly to ID arrays we find, that one 
impurity is able to organize the 2D arrays. However, our 
bifurcation analysis shows that the coupling constant has 
to be larger than fccr — 3, whereas in the ID case k > 0.1 
was sufficient fl^ to produce spatio-temporal patterns. 
We would like to point out the appearance of a P4 at- 
tractor for k — 5, Innp < 0.15 in Fig. |l]i. Here an array 
of chaotic pendula is synchronized by a single impurity, 
which itself is chaotic, when isolated. 
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FIG. 1. Bifurcation diagram for a 50 x 50 lattice of coupled 
pendula with length Zn,m = 1- A single impurity with length 
limp is located at lattice site (25,25). For each hmp the values 
of ct(151T), ...,(t(170T) are shown, where T is the period of 
the driving. The coupling is (a) k — 1, (b) k = 2, (c) fe = 3, 
and (d) fc = 5. In (c) the first windows of synchronization 
can be observed, which enlarge for bigger couplings (d). 



In order to illustrate the occurence of this P4 pattern 
in a better way, we show in Fig. ^ a typical gray scale 
plot of the velocities as a function of the lattice 

coordinates (n, to). The stroboscopic snapshots are taken 
with a time lag of one period of the driving T; darker 
shading indicates higher velocity. In this example we 
used 64 X 64 oscillators having Z„^„i — 1 and located an 

0.1 in the middle of the lattice 
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impurity of length I 
at (32,32). The coupling constant is A; = 5. After some 
transient (not shown) the synchronization is maintained 
and a P4 attractor is clearly visible. We neglect edge 
effects except to remark that they are strictly confined 
to the last few pendula near the boundaries. Moreover, 
they become very regular if the calculation is carried on 
for longer times. 
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FIG. 2. Snapshots of 6„^m for a 64 x 64 lattice with 



1, hn 



imp — ^32, 32 

of the driving period T. 



0.1 and A: = 5 at multiple times 



In many applications, however, one is interested in the 
weak coupling limit p6|. In this limit, our analysis shows 
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that one impurity is not able to create spatio-temporal 
organization of a chaotic lattice. Therefore, it is meaning- 
ful to ask whether spatio-temporal patterns can emerge 
at all and under which conditions they may be obtained. 
To this end we increase the number of impurities intro- 
duced in the lattice and investigate the importance of 
their geometrical arrangement as well as their concentra- 
tion. 

In Fig. |3^ we report the PI pattern emerging for a 
128 X 128 lattice with coupling constant k = 0.5, when 
128 impurities were arranged along a line in the middle 
of the lattice, dividing it into two disjoint parts. The 
same behavior could be observed even for smaller cou- 
plings. We verified this for coupling constants as small 
as fccr — 0.1. Moreover, we found that increasing the 
size of the lattice to TV = 256 (the limit of our computa- 
tional capability) , but maintaining the line- like geometry 
of impurities and the parameter-values, did not affect the 
formation of a pattern. In all cases we observed synchro- 
nization to a PI pattern after an initial transient. Thus, 
increasing the size of the lattice will not affect the pat- 
tern formation although the percentage of the impurities 
will tend to zero in the limit of infinite systems. This 
indicates that the concentration of impurities is not of 
primary importance. 

To test the influence of the special arrangement of im- 
purities, we considered a 128 x 128 lattice with k = 0.5 
and 128 impurities of length Zi,„p — 0.7 located at ran- 
dom positions of the lattice. In all cases we have studied 
we obtained chaotic patterns. Even an increase of their 
concentration by more than a factor of 3 did not produce 
any spatio-temporal pattern. 
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FIG. 3. Snapshots of ^„,m for a 128 x 128 lattice with 
ln,m ~ 1, 'imp = 0.7 and k = 0.5. The impurities are po- 
sitioned as indicated in the upper panel: (a) on a line along 
n = 64, (b) on a cross along the lines n — 64 and m — 64, (c) 
on a square with (n, m) = (63, 63) the lower left and (66, 66) 
the upper right corner. In the lower panel, the corresponding 
(y{jT) is shown. 



It is thus natural to ask wether a line of impurities 
is the only geometry that can produce spatio-temporal 
patterns for relatively small values of the coupling con- 
stant. In Fig. ||b,c we report the PI attractors of a lat- 
tice of 128 X 128 coupled pendula with lengths ln,m = 1, 
hmp = 0.7 and coupling k = 0.5, where the impurities are 
arranged like a cross (Fig. ^d) or as a ring (Fig. He). In 
both cases the observed pattern geometries consist locally 
of stripes as was also observed for the 'line-geometry' in 
Fig. |a. 

What is the common geometrical feature of the above 
impurity configurations that allow them to control the 
chaotic lattice? From the above analysis we draw the 
conclusion, that it is sufficient for synchronization, that 
the impurites divide the lattice into at least two disjoint 
subregions. In that way the critical coupling needed to 
observe a spatio-temporal pattern is decreased by an or- 
der of magnitude. Reducing the coupling below fccr ~ 0.1, 
however, suppresses the formation of spatio-temporal 
patterns even for these lattices. This finding is consistent 
with earlier investigations on ID arrays, which revealed 
a critical coupling k « 0.1, below which no pattern for- 
mation could be observed. Furthermore, the observed PI 
pattern for the 'line geometry' (Fig. ^) is analogous to 
the ID case, where the introduction of a single impurity 
caused the same topology of the array, i.e. a division into 
two disoint sets, and a similar pattern was observed p4[ . 
Thus, the results of the ID case define the limiting fc- 
value also for the 2D lattice. 
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FIG. 4. (a) Chaotic sequence a{jT) for a 128 x 128 lat- 
tice (k = 0.5) with e [0.998, 1.002]. (b) Periodic 
sequence o-(jT) for a 128 x 128 lattice (k = 0.5) with 
ln,m € [0.8, 1.2]\[0.998, 1.002] revealing a P6 pattern, (c) 
Snapshots at multiple periods T of the driving confirm the 
existence of a P6 pattern as predicted in (b). 



We finally examined the effect of disorder and the pos- 
sibility of obtaining self-organization or frequency lock- 
ing. For a 128 x 128 lattice we randomly varied the 
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lengths of the pendula but restricted the range of the 
disorder such, that each individual pendulum is chaotic, 
i.e. ln,m e [0.998,1.002]. We found that the emerg- 
ing pattern was always chaotic for a large number of 
different realizations of disorder and initial conditions. 
The lacking synchronization can be observed in a bet- 
ter way by inspecting cr(jT), which does not show any 
periodicity (see Fig. ^ for an example). If the range 
of disorder is increased to include also regular moving 
pendula, i.e. ln,m G [0.8,1.2], self organization is possi- 
ble in agreement with findings of Ref . jl2| . This can be 
understood by the observation that synchronization al- 
ready occurs, when choosing the disorder from an inter- 
val of pendulum lengths associated with regular motion, 
i.e. ln,m e [0.8, 1.2]\[0.998, 1.002]. A particular example 
for this is shown in Figs. ^,c, where the occurence of 
a P6 pattern can be observed. Thus, taking the lengths 
from the entire interval ln,m, G [0.8, 1.2] on average yields 
only 1% chaotic pendula, whose motion will be overdom- 
inated by the synchronizing regular ones, explaining the 
spatio-temporal pattern observed in Ref. |l2| . 

IV. CONCLUSION 

In conclusion we have demonstrated that a lattice of 
chaotic pendula can be frequency locked into a spatio- 
temporal pattern by introducing impurities in the lattice. 
In the strong coupling limit, a single impurity can tame 
chaos. Decreasing the coupling constant requires more 
impurities in order to observe self-organization. In this 
case the geometry of the impurity configuration plays an 
important role. Our results suggest that if the impurity 
configuration divides the lattice into at least two disjoint 
parts then the coupling constant may be decreased with- 
out affecting the synchronization of the chaotic array. 
Moreover, the induced spatio-temporal patterns are then 
unaffected by the size of the lattice. Below the critical 
coupling fccr ~ 0.1 no synchronization is observed. The 
value of this critical coupling is dictated by the minimum 
coupling which leads to the formation of spatio-temporal 
patterns in the ID case. 

We are grateful to T. Gavrielides, V. Kovanis, A. Politi, 
and G Tsironis for helpful discussions. 
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